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We calculate the screening function in bilayer graphene (BLG) both in the intrinsic (undoped) and 
the extrinsic (doped) regime within random phase approximation, comparing our results with the 
corresponding single layer graphene (SLG) and the regular two dimensional electron gas (2DEG). 
We find that the Kohn anomaly is strongly enhanced in BLG. We also discuss the Friedel oscillation 
and the RKKY interaction, which are associated with the non-analytic behavior of the screening 
function at g = 2kF- We find that the Kohn anomaly, the Friedel oscillation, and the RKKY 
interaction are all qualitatively different in the BLG compared with the SLG and the 2DEG. 
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Single layer graphene (SLG), a single layer of carbon 
atoms arranged in a honeycomb lattice, has attracted a 
great deal of attention, both experimentally and theoret- 
ically, for its unusual electronic transport and the charac- 
teristics of relativistic charge carriers behaving like mass- 
less chiral Dirac fermions [1]. Bilayer graphene (BLG) 
consisting of two SLG is also of great current interest, 
both for technological applications and fundamental in- 
terest [l|, [2, _3J . While the band structure of SLG has a 
linear dispersion, BLG has a quadratic dispersion ^2| in 
the low energy regime making it similar to two dimen- 
sional (2D) semiconductor systems except for the absence 
of a gap. The purpose of this work is to calculate the 
polarizability (or screening) function of bilayer graphene 
within the random phase approximation (RPA). Even 
though many theoretical works on Coulomb screening in 
SLG have been reported [J, [sj, ISJ, [Zl, |8| , the analytic in- 
vestigation of Coulomb screening in BLG has not yet 
been performed. Knowing the BLG screening function is 
crucial since it determines many fundamental properties, 
e.g. transport through screened Coulomb scattering by 
char ged impurities [9| , Kohn anomaly in phonon disper- 
sion [lO|, and RKKY interaction 11|. In order to un- 
derstand the electronic properties of BLG it is therefore 
necessary to obtain its screening function. 

The BLG is in some sense intermediate between the 
SLG and the regular semiconductor-based two dimen- 
sional electron gas (2DEG) since it is chiral with a zero 
band gap at the Dirac point (where the electron and 
the hole bands touch) similar to the SLG, but has the 
quadratic energy dispersion similar to the 2DEG. For ex- 
ample, the 2fci?-backscattering is suppressed y, la, la, [Zl, l8| 
in SLG due to its chiral nature whereas in the 2DEG, the 
2A:F-backscattering plays a key role 112] in determining 
low density and low temperature carrier transport. We 
find that the 2A:F-backscattering is restored (and even en- 
hanced) in the BLG because of the quadratic dispersion 
and, more importantly, due to the symmetry imposed by 
the two-layer structure. This qualitative difference in the 



screening properties between BLG and SLG leads us to 
predict that transport and other electronic properties in 
BLG will be qualitatively more similar to 2DEG than to 
SLG in spite of the zero-gap chiral nature of BLG. 

The effective BLG Hamiltonian is now well established 
in the theoretical literature. In the low energy regime the 
Hamiltonian is reduced to the (2 x 2) matrix form and is 
given by (we use h = 1 throughout this paper) [2] 
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where m = ji/{2vp), 71 is the interlayer tunneling anrpli- 
tude, and vp is the SLG Fermi velocity. The wave func- 
tion of Eq. H]) can be written as ipsk = e*''(e~2^^'', s)/V2 
and the corresponding energy is given by esk = sk'^ /2m, 
where 0k — ta.n^^ (ky/kx) and s = ±1 denote the band 
index. Using the Hamiltonian of Eq. ([T|) we theoretically 
obtain the screening function of BLG by calculating the 
polarizability and the dielectric function within RPA. 
The static dielectric function can be written as 
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where k is the background dielectric constant, and H(g) 
the polarizability. The static BLG polarizability is given 
by the bare bubble diagram 
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where g is the degeneracy factor (here .g = 4 due to valley 
and spin degeneracies), k' = k + q, s,s' = ±1 denote 
the band indices, Esk = sk'^/2m,, and i^ss'(k, k') = (1 -|- 
ss' cos26')/2, where 9 is the angle between k and k', and 
fsk is the Fermi distribution function, /^k = [exp{/3(esk — 
/i)} + l]~^, with f3 = l/ksT and /i the chemical potential. 
First, we consider intrinsic (i.e. undoped or ungated, 
with n and Ep both being zero) BLG where the conduc- 
tion band is empty and the valence band fully occupied 



at T = 0. Then we have /k+ — and /k- = 1. Since 
the conduction band is empty the polarization is induced 
by the virtual interband transition of electrons from the 
valence to the conduction band. The polarizability due 
to the interband transition becomes n(g) = n''((7), where 
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where cos6' = (|k| 
calculated easily 



£+k — E-k' E-k — e+k' 
q|cos0)/|k + q|. Eq. (g]) can be 
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(5) 



where iVp — gm/2n is the BLG density of states. Thus 
the intrinsic BLG polarizability is constant for all wave 
vectors. (Note that the polarizability of ordinary 2DEG 
is constant [11[ only for q < 2kp.) The dielectric function 
becomes e(g) — l + qg/q., where the screening wave vector 
is given by 



qs = qTF log 4, 



(6) 



whe re q tf is the 2D Thomas-Fermi screening wave vec- 
ll| . qTF = grne^JK. BLG static screening is thus 



tor 



enhanced by a factor of log 4 compared with ordinary 
2D screening till]. For intrinsic BLG we can write the 
screened Coulomb potential as 
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where Hq(x) and Nq(x) are the Struve function and the 
Bessel function of the second kind, respectively. The 
asymptotic form at large r is 4'{r) ^ eqs/K{qsr)^- Since 
the screening function is a constant for all q without any 
singular behavior there is no oscillatory term in the po- 
tential. This is very different from the screening behavior 
of intrinsic SLG [ISli] or 2DEG [uj. 



For intrinsic SLG we have 
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where N^^'^ = gq/{2TrvF)- The intrinsic SLG po- 
larizability increases linearly with q, and e{q) = 1-1- 
{e'^g/KVp){Tr/8), which gives rise to only an enhance- 
ment of the effective background dielectric constant k* — 
K + {e'^g/vp){-K/%), i.e. the screened Goulomb interac- 
tion V{q) — 2-^6^ / K*q. The Coulomb interaction in real 
space can be expressed by V{r) — e^ / K*r for all r. Thus 
at large r the Coulomb potential decreases as 1/r^ in 
intrinsic BLG, but only as 1/r in intrinsic SLG. 

In the following we provide the zero temperature static 
polarizability of extrinsic (i.e. gated or doped) BLG 
where n,EF ^ 0. At T=0, /_k = 1 and /+k = 
0{kp — |k|). Then we can rewrite Eq. (1) as n((7) = 

nintra(9) + niiitor(g), whcrC 
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Hintra (Hintcr) indicates the polarization due to intraband 
(interband) transition. After angular integration over the 
direction of q we have 
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Finally, we have the extrinsic BLG static polarizability 
as 



n[q)^No[J{q)~g{q)e{q~2kp)]. 



(16) 



Eq.(16) with Eq. (15) is the basic result obtained in this 
paper, giving the doped BLG polarizability analytically. 
In Fig. [1] we show the calculated static polarizabil- 
ity as a function of wave vector. Fig. [1] (a) and (b) 
show the calculated intraband and interband polarizabil- 
ities, respectively, with those of single layer graphene for 
comparison. Fig. [1] (c) shows total polarizability of bi- 
layer graphene. At g = we have nintra(O) = A^o and 
nintor(O) = which follow also from the compressibil- 
ity sum rule, n(g = 0) = J de {-df{e)/de) N{eg). For 
small q, nintra(g) dccrcases as 1 — q'^/2kp, and nintcr('z) 
increases as q'^ /2kp. This behavior comes from the over- 
lap factor Fss' in Eq. ^. For SLG intraband (inter- 
band) polarizability decreases (increases) linearly as q 
increases and these two effects exactly cancel out up to 



q = 2fci?, which gives rise to the total static polarizabiUty 
being constant for q < 2kF as in the 2DEG. However, for 
BLG the cancehation of two polarizabihty functions is 
not exact especiahy for q> kp because of the enhanced 
backscattering, so the total polarizabihty increases as q 
approaches 2k p, which means screening increases as q 
increases. Thus BLG, in spite of being a 2D systems, 
does not have a constant Thomas-Fermi screening up to 
q — 2kp as exists in SLG and 2DEG. 

A qualitative difference between SLG and BLG po- 
larizabihty functions is at g = 2kF- Due to the sup- 
pression of 2kF backward scattering in SLG, the total 
polarizabihty as well as its first derivative are contin- 
uous. In BLG, however, the large angle scattering is 
enhanced due to chirality, which gives rise to the singu- 
lar behavior of polarizabihty at q — 2kp. Even though 
the BLG polarizabihty is continuous at q = 2kF, it 
has a sharp cusp and its derivative is discontinuous at 
2k p, diverging as q approaches 2k p, i.e., as q —^ 2k p, 
dli{q)/dq (x l/y^q"^ — 4:kp. This behavior is exactly the 
same as that of the regular 2DEG, which also has a 
cusp at g = 2kp in addition to being constant in the 
< g < 2kp region. Note that in SLG this non-analytic 
behavior of polarizabihty occurs in the second derivative, 
d^U{q)/dq^ OC l/y'q^-Akj,. 

In the large momentum transfer regime, q > 2kp, the 
BLG polarizabihty approaches a constant value (intrin- 
sic polarizabihty 11°), i.e., n(g) -^ A'^olog4, because the 
interband transition dominates over the intraband con- 
tribution in the large wave vector limit. This is very 
different from that of a 2DEG where the static pa 



larizability falls off rapidly (~ 1/g^) for q > 2fci?[ll| 



and SLG where the polarizabihty increases linearly as q. 
Thus, the large q behavior of dielectric screening becomes 
e(g —> oo) — > l+gTre'^/{8KVp) for SLG and e{q -^ cx)) — > 1 
for both BLG and 2DEG. 

The strong cusp in BLG n(g) at g = 2kp leads to 
Friedel oscillations in contrast to the SLG behavior. The 
leading oscillation term in the screened potential at large 
distances from a point charge Ze can be calculated as 
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where C = \/5 — log[(l + ^/b)/2], which is similar to the 
2DEG except for the additional constant C (C = 1 for 
2DEG), but different from SLG where Friedel oscillations 
scale as (/)(r) ~ cos(2fci?r)/r'^. The enhanced singular 
behavior of the BLG screening function at g = 2kp has 
other interesting consequences related to Kohn anomaly 
|ld | and RKKY interaction, which we discuss below. 

The strong cusp in the BLG polarizabihty at g = 2k p^ 
as can be clearly seen in Fig. 1(c), indicates that the 
screened BLG acoustic phonon frequency would mani- 
fest a strong Kohn anomaly, i.e. an observable dip struc- 
ture in the phonon frequency at g = 2kp. It is obvious 
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FIG. 1: (Color online) Calculated (a) intraband, (b) inter- 
band, and (c) total static polarizability of bilayer graphene. 
For comparison the single layer polarizabilities are shown. In 
(c) we also show the regular 2D static polarizability (dashed 
line) . 



from Fig. 1(c), and from the discussion above based on 
our analytical results for e(g), that the screened phonon 
dispersion will exhibit a much stronger Kohn anomaly 
in the BLG than in the SLG with the 2DEG being in- 
termediate. This arises from the stronger singularity at 
g ~ 2kp manifesting in the first derivative dll{q)/dq in 
the BLG rather than in the second derivative, d'^H/dq^, 
in the SLG. Tuning the value of kp by changing the car- 
rier density through the applied gate voltage, it should 



be possible to verify that the Kohn anomaly is indeed as- 
sociated with the 2kF screening behavior in the BLG. In 
fact, the BLG Kohn anomaly would be, due to the very 
strong q = 2kp cusp in the polarizability, rather simi- 
lar to the ID Peierls instability [13| since the q — 2kp 
screening behavior in BLG is qualitatively similar to the 
ID electron system '14']. 

The polarizability function in Eq. (fT6|) also deter- 
mines the RKKY interaction between two magnetic im- 
purities due to the induced spin density (here we con- 
sider magnetic impurities located at the interface be- 
tween BLG and substrate, so they do not break any sym- 
metry). The RKKY interaction (or induced spin den- 
sity) is proportional to the Fourier transform of ^{q)- 
The conventional form of the exchange interaction be- 
tween the localized moment S and electron-spin density 
s(r) is given by V{r) — JS'(R)s(r)5(R — r), where J 
is the exchange coupling constant. The RKKY inter- 
action between two localized moments via the conduc- 
tion electrons may then be written in the following form: 
HRKKvir) = J^5i5'2n(r), where U{r) is the Fourier 
transform of the static polarizability n(g). 

First consider intrinsic BLG. The Fourier transform of 
Eq. (O simply becomes a (5- function because Ii°{q) is a 
constant, i.e. n(r) = NQ\og{4)5{r). This indicates that 
the localized magnetic moments are not correlated by 
the long range interaction and there is no induced spin 
density. In SLG, the Fourier transform of polarizability 
(Eq. ([8])) diverges (even though n(r) formally scales as 
1/r^, its magnitude does not converge), which means that 
intrinsic SLG is susceptible to ferromagnetic ordering in 
the presence of magnetic impurities |15l . Il6l | due to the 
divergent RKKY coupling. 

In doped (or gated) BLG the oscillatory term in RKKY 
interaction is restored due to the singularity of polariz- 
ability at 5 = 2kp, and the oscillating behavior domi- 
nates at large kpr. At large distances 2kpr ;» 1, the 
dominant oscillating term in 11 (r) is given by 



n(r) - iV( 



kp sm{2kpr) 
2^ {kprf 



This is the same RKKY interaction as in a regular 2DEG, 
and it decreases as l/r^, in contrast with 1/r^ behavior 
in SLG. @ 

In conclusion, we calculate analytically the static wave 
vector dependent polarizability of both undoped and 
doped bilayer graphene within RPA. For undoped BLG 
we find that screening is enhanced by a factor of log 4 
compared with ordinary 2D screening. The RKKY in- 
teraction in undoped BLG is zero-ranged (^-function) 
and therefore, no spin density is induced. The doped 
BLG screening function shows strongly enhanced Kohn 
anomaly at 2kF compared with the corresponding SLG 
and 2DEG situations, which give rise to the usual RKKY 
interaction and Friedel oscillation. We show that BLG 
screening properties are qualitatively different from SLG 



screening behavior in all wave vector regimes {q < 2/cf, 
q > 2kp, and q = 2kp) with the BLG screening having a 
strong cusp a,t q = 2kp. Our theory applies only in the 
low density regime (n < 5 x lO^^cm"^), where the band 
dispersion is quadratic and only the lowest subband is oc- 
cupied [4] . There are obvious implications of our results 
for BLG carrier transport limited by screened Coulomb 
scattering - in particular, the strong 2kp anomaly in 
screening will lead to strong temperature dependence in 
dc transport at low (T << Tp) temperatures. This is 
in sharp contrast to SLG where 2kp backscattering is 
suppressed. 
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